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Abstract 

We present a proof of concept design for a structured square cloak enclosing a void in 
an elastic lattice, subjected to flexural vibrations. In addition to the theoretical design, 
we fabricate and experimentally test an elastic invisibility cloak for flexural waves in a 
mechanical lattice. Eurther verification and modelling is performed numerically through 
finite element simulations. The primary advantage of our square lattice cloak, over other 
designs, is the straightforward implementation and the ease of construction. The elastic 
lattice cloak shows an impressive efficiency within the predicted frequency range of 100- 
250Hz. 
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1 Introduction 

Since the first experimental demonstration of the microwave invisibility cloak [1], there has been 
an explosion of theoretical and practical advances in the design and analysis of electromagnetic 
metamaterials [2]. In contrast, the much more challenging problem of creating invisibility cloaks 
and metamaterials for elastodynamics has been much less studied. Notable recent advances in 
the theoretical analysis of cloaks for elastic waves have been made by Milton et al. iiii], Norris 
&; Shuvalov [5], Brun et al. lao, Farhat et al. [8], Jones et al. [9], Colquitt et al. uniin], 
Guennau et al. [12] and Parnell et al. [nmansiiiB]. These theoretical developments have been 
complemented with a series of experimental implementations of multi-scale mechanical cloaks 
performed by the group led by Wegener na ng [H]. 

In particular, the work by Milton, Briane, &: Willis 1311 ] identifies the effects of ‘negative 
inertia’ in the elastic cloak and analyses the requisite strong anisotropy, not only in the elastic 
compliance, but also in the inertial properties of the cloak. The fact that, in the framework of 
Milton et al. 1313 , the mass density of the material should behave as a tensorial quantity rather 
than simply as a scalar is an important and striking observation. A novel approach of dynamic 
homogenisation, capable of encapsulating such striking effects, was introduced and systemati¬ 
cally studied by Craster and co-authors mi EH- The asymptotic theory leads to an effective 
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Fig. 1: Two-dimensional ABAQUS simulations performed in the case of a plate without a hole (upper), a plate 
with an uncloaked hole (middle), and a plate with a hole surrounded by the cloak (lower). For a frequency range 
between 120 Hz and 240 Hz, the simulations show that the cloak reduces significantly the scattered field. 


equation for the envelope function in the perturbation approximation relative to a standing 
wave associated with a periodic lattice or continuum structure. Norris & Shuvalov [5] later 
generalised the framework of Milton, Briane, & Willis and showed that, although one cannot 
create an invisibility cloak for elastodynamics without recourse to a non-classical generalised 
theory of elasticity, with an appropriate choice of Gauge one can choose between creating a mi- 
cropolar elastic cloak or a cloak with tensorial density, for example. The concept, design, and 
theoretical analysis of elastodynamic invisibility cloaks are much more challenging compared to 
cloaks for membrane, acoustic, optical, and anti-plane shear waves, all of which are governed by 
the transformed Helmholtz equation. Although the theoretical framework for elastic cloaks is 
well defined, its experimental implementation has never been successfully achieved for dynamic 
vector problems of elasticity. The pioneering work by Wegener and his group for Kirchhoff-Love 
plates nzi is the only significant experimental contribution in this extremely challenging area. 
In a different context, multi-scale resonators were discussed in j22] in relation to an approximate 
cloaking referred to as ‘Directional cloaking’ for elastic plates containing voids. 

The elastic Kirchhoff-Love plate provides an efficient and rigorous framework for the analysis 
of elastic waves in thin plates. The propagation of flexural waves in a Kirchhoff-Love plate is 
governed by the biharmonic operator and waves in an homogeneous isotropic plate can be 
expressed as a linear combination of solutions of the Helmholtz equation (which we refer to as 
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Fig. 2: Comparison between experiments and numerical simulations in the case of an applied displacement with 
a frequency of 120 Hz 


“membrane waves”) and solutions of the modified Helmholtz equation. The latter are evanescent 
fields but, nevertheless, may make a significant contribution through the boundary conditions 
and hence play a crucial role in the dispersion of flexural waves in structured plates [saisiES], 
in addition to the design of flexural invisibility cloaks [9]. 

The design of the cloak presented here is distinct from the cloaks presented by Wegener 
et al. m, in so far as the cloak designed and constructed here is matched with a multi-scale 
mechanical structure, as opposed to a homogeneous continuous ambient matrix. Moreover, the 
cloak which we construct here is not singular, i.e. it does not require infinite wave speeds on the 
interior boundary, but is based on the regularised cloaking transformation and theoretical design 
presented by Colquitt et al uni for membrane waves and later for flexural waves in plates m- 
The cloak design is shown in Fig. and the numerical and experimental implementations 
demonstrate that the invisibility cloak is efficient. This is illustrated in Fig. [T| which shows a 
uniform vibrating plate without a hole, a plate with an uncloaked hole, and a plate with a hole 
surrounded by the cloak - the latter reduces significantly the scattered field, as predicted by 
the analysis. 

It is emphasised that here we consider the full elastodynamical problem of wave propagation 
in a discrete metamaterial lattice cloak embedded within a multi-scale ambient medium; this 
should be distinguished from the earlier work of Wegener et al. on mechanical lattice cloaks m 
uni wherein the static cloaking problem is considered and no waves propagate within the system. 
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The structure of the paper is as follows. Sectionj^outlines the idea of the regularised cloaking 
transformation leading to the practical implementation of the multi-scale structured square 
cloak, which is surrounded by an ambient square lattice. The computational and experimental 
implementations are discussed in Section]^ In particular, the computational section |TT] presents 
several ABAQUS simulations to verify the efficacy of the cloaking effect, as shown in Fig. 
The experimental design and results are discussed in section |4.2[ Finally, we draw together 
important concluding remarks in section 

2 The Hooke-Chladni-Faraday visualisation 

Before presenting the mathematical model of the cloak and the experimental implementation, 
we first discuss an elegant technique that we will employ to visualise the time-harmonic vibration 
of plates. The method, known as the ‘^Hooke-Chladni-Faraday’^ technique, has been employed 
by many researchers in physics and mechanics for almost four centuries and continues to be 
used today. Indeed, the recent paper [26] examines the visualisation of standing waves in 
dynamically reconfigurable liquid-based templates, which were used to assemble micro-scale 
materials into ordered structures with desired properties. This novel approach at a micro-scale 
level employed the idea of Faraday waves. A three-dimensional visualisation of acoustic standing 
waves was reported in a recent paper m that demonstrated the elegant efficiency of the classical 
Hooke-Chladni-Faraday method, which continues to generate new ideas and exciting results. In 
particular, the experiment reported in m, shows levitating micro-particles along the nodal 
lines of a three-dimensional standing wave. 

Well before the time of Michael Faraday, Robert Hooke and then Ernst Chladni discovered 
an ingenious method to visualise standing waves in elastodynamics. This was especially effective 
for flexural resonances in elastic plates and membranes, as described in Chladni’s book [28]. The 
technique consists of drawing a violin bow over a metallic or glass plate that is covered with 
flour. Once the plate reaches a resonance the flour collects along the nodal lines of the resonant 
mode providing an elegant yet efficient visualisation of the standing wave present in the plate. 

Following Chladni’s experimental demonstration, the eigenvalue problem for the free vibra¬ 
tions of a square plate with free edges has been studied by many scientists, most notably. Lord 
Rayleigh [23 ED] and Ritz 1311E2] during the development of the now well known Rayleigh- 
Ritz method. Although no closed form solution currently exists, Ritz [32], over a century ago, 
was able to construct remarkably accurate approximate solutions for the eigenfrequencies and 
nodal patterns. These classical studies have led to a fascinating area of Cymatics, which anal¬ 
yses methods of making sound and vibration visible and has attracted attention of engineers, 
mathematicians, physicists and musicians around the world. 



Fig. 3: Hooke-Chladni-Faraday visualization of four eigenmodes of a square elastic plate with a free boundary. 
These illustrative experiments were produced in the ‘Instabilities Lab’ of the University of Trento. 
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For the purpose of illustration, we include in Fig. examples of the Chladni patterns 
for eigenmodes of a square elastic plate with a free boundary, which are accurate and fully 
consistent with analytical findings. The patterns depend on the boundary conditions and on 
any inhomogeneities that may occur in the plate, such as voids or inclusions. In particular, 
Chladni patterns were never constructed for a plate with a hole surrounded by a structured 
cloak. 

In the present paper, we visualise Hooke-Chladni-Faraday patterns for flexural waves around 
an obstacle surrounded by a multi-scale structured cloak and elegantly illustrate the efficacy of 
the cloak. An ABAQUS simulation for the mechanical configuration, identical to the one used 
in the experiment, provides accurate numerical data on the displacement amplitudes and stress 
distribution. We present the experimental visualisation to confirm the predicted wavefront 
profiles. 

Compared to the classical settings for the above mentioned frequency response problem for 
a rectangular Kirchhoff plate, we go further and consider a structured plate with a cloaked 
hole. Fig. includes Hooke-Chladni-Faraday patterns used for visualisation of the cloaking 
effect of flexural waves. These observations are new and demonstrate scattering patterns for 
three configurations, which include (a) a rectangular lattice-type plate, (b) a plate with a square 
hole, and (c) a plate with a structured cloak enclosing the hole. As in the original experiments 
by Hooke, the powder collects along the nodal lines of the vibrating plate thus indicating the 
position of the wavefronts, i.e. the locus of points on the wave with the same phase and 
zero displacement. The Hooke-Chladni-Faraday patterns allow us to conveniently visualise the 
propagation of waves in the structured plate and, in particular, we can observe a significant 
reduction of the scattering pattern for the case (c) (as in the bottom part of Fig. [^. A detailed 
discussion of the experiment is given in the main text of the paper. 

3 The square cloak 

In most of the papers addressing the design of cloaks for linear waves, a singular radially sym¬ 
metric push-out transformation is employed (see, for example, 1331 USD which maps a point to 
a finite disc. For theoretical and computational models of continuous media, such an approach 
is adequate. However, for practical implementations it poses substantial difficulties. In partic¬ 
ular, not only does the transformation lead to infinite wave speeds on the interior boundary 
of the cloak, but also the required cloaking material is characterised by an unrealistic strong 
anisotropy. Indeed a lattice structure, instead of a continuum, can be used to create an in¬ 
visibility cloak. A significant advantage of lattices is that they naturally accommodate high 
contrasts in their compliance leading to strong anisotropy. However, a radial discrete cloak that 
fits inside a circular ring would not match any periodic lattice, and the presence of a geometrical 
mis-match on the interface boundary leads to a substantial mis-match in the interface boundary 
conditions. 

For the design presented here, we choose to employ a square cloak following the theoretical 
framework recently established in m- The cloak is embedded in an ambient square lattice, 
which is subjected to out-of-plane flexural vibrations. We consider a thin structured plate with 
a defect represented by a traction free void, which significantly influences the wave field, as 
illustrated in Fig. A structured cloak, as in [iniiii], is then installed around the void. In the 
cloaked configuration, we observe a significant reduction in the scattered field and, in particular, 
the reduction in the shadow region behind the scatterer and the restoration of the incident field 
represented by a plane wave. 
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Fig. 4: (A) The function maps the undeformed trapezoidal region (1) to the deformed trapezoidal conhgu- 
ration. (B) A discrete lattice structure where the curved ligaments are oriented following the principal directions 
of the stiffness matrix for the continuum cloak. 


3.1 The regularised cloaking transformation 

Following Colquitt et al. we choose a regularised near-cloak, which is obtained by four 

push-out transformations applied to trapezoidal regions, as illustrated in Fig. ih- The idea 
of regularisation is to map a domain with a small hole (e.g. a square of semi-width e as in 
Fig. at the centre into another domain, whose exterior boundary is preserved while the 
interior boundary is expanded to the required finite size. For the regularised problem, we set 
the Neumann boundary condition on the boundary of the hole, which means (in the case of an 
elastic plate) the free edge boundary condition. As showed recently [9], the correct choice of 
boundary condition is of vital importance in order to achieve cloaking. 

The mapping is defined in such a way that x = for each trapezoidal region {j = 

1,2, 3,4) shown in Fig|^. In particular, the map for the trapezoidal region (1), with bases 
perpendicular to the Xi—axis, is defined by the formula: 




CXlXi -|- CX2 

0^1X2 + a 2 X 2 /Xi 


where 


w 


OL\ — 


0L2 — 


{a + w) {a — e) 


( 1 ) 


( 2 ) 


a -\- w — e a -\- w — e 

Here, 0 < e 1 is the regularisation parameter and uppercase letters denote the undeformed 
configuration, whilst lowercase letters denote the deformed configuration. Prior to the transfor¬ 
mation, the interior boundary of the trapezoidal region (1) corresponds to Xi = e. Following 
the transformation, the interior boundary of this trapezoidal region is moved to xi — a. The 
outer boundary of the cloak is invariant with respect to the transformation, and for the region 
(1), it corresponds to Xi = a + w. The other three trapezoidal regions are transformed in the 
same way, subject to a rotation. 


3.2 Transformed equations after the cloaking transformation 

The time-harmonic flexural deformations of a homogeneous thin elastic plate are governed by 
the Kirchhoff plate equation, which before the transformation has the form 

{v%-^u^)w{X) = 0. 
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After the application of a non-singular transformation we have 

(^V • J-^FF^VJV • J-^FF^V - ^Lo‘^^w{x) = 0, (3) 

where F = is the non-degenerate Jacobi matrix of the geometrical transformation, and J — 
det F. The above equation can be interpreted as the equation of an anisotropic pre-stressed plate, 
as discussed in m- In the context of asymptotic approximations, the paper [7] discusses the 
configurations where the effects of pre-stress are small, and when an approximate factorisation 
of the transformed differential operator is possible. 

Our goal is to construct an approximate cloak, which may be conveniently implemented 
experimentaly. We observe that in the successful experiment for the continuous plate led by 
Wegener nzi, the authors assumed that the membrane waves are dominant and configured 
their radially symmetric approximate cloak accordingly. Although such an assumption may 
appear to be inappropriate, it has been demonstrated that, within a certain frequency range, 
the membrane waves indeed are dominant and the approximate cloaking effect is apparent. 
This approach is further reinforced for the square cloak by the observation made in im that 
the principal directions of stiffness for the membrane cloak and the flexural cloak are exactly 
the same. Compared to implementations of invisibility cloaks for continua, such as in 
the square lattice cloak developed in the present paper is relatively straightforward. We now 
proceed to describe the implementation of our novel square lattice cloak. 

3.3 The lattice approximation of the cloak 

We consider an approximate cloak, realised using a discrete lattice structure with curved liga¬ 
ments as derived in m- These elastic ligaments are aligned with the the principal directions 
of the stiffness matrix for the continuum cloak, as illustrated in Fig. P, and yield a structured 
flexural lattice system. The thin elastic ligaments may be treated as beams of rectangular 
cross section characterised by the bending stiffnesses Di and D 2 chosen in accordance with the 
analytical formulae from im 


Di 


Xi — OL2 

5 

XX 


xf -h Oi^x\ 
x\{xx — 0 ^ 2 ) 


(4) 


These flexural rigidities are visualised in Fig. [^for the right-hand quadrant. For this quadrant, 
on the Xx axis, Dx represents the Xx stiffness in the principal direction, whereas D 2 represents 
the X 2 stiffness in the other principal direction. In every other point inside the cloak, Dx and D 2 
represent the stiffnesses in the principal directions of the locally orthotropic cloak constructed 
here. We note that on the interior boundary of the cloak the tangential rigidity D 2 is much 
higher compared to the normal rigidity Dx and emphasise that these stiffnesses are finite on 
this boundary of the cloak (see Fig. [^. 

According to equation i> the mass density p' inside the cloak is also non-uniform and obeys 
the formula 


P 


/ 


Xx ~ 0^2 

P 2 ’ 

a^xx 


(5) 


where p is the density outside the cloak and the coefficients ax and 0^2 are defined in Q. 

Equations @ and together with the accompanying geometrical design of the cloak, 
provide the essential information for the experimental implementation discussed in the text 
below. 
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Fig. 5: The required stiffnesses Di and D 2 for the cloak ligaments reported as function of xi and X 2 . 


4 Numerical simulation and experimental verification 

We have performed both finite element (FE) computations, in ABAQUS, and also experiments 
to verify the efficacy of the square structured cloak in reducing the scattering of flexural waves 
from voids. The structured cloak has been designed and implemented in SOLIDWORKS. Each 
elastic ligament has a specified variable cross section that provides the required rigidities Di 
and D 2 as reported in Colquitt et al m- The physical lattice cloak was created by milling 
holes into a polycarbonate plate; both ABAQUS and the milling machine were programmed 
using the same SOLIDWORKS original code. In so doing, experiments and simulations have 
identical geometry, material parameters, constraints and applied out-of-plane vibrations. Figure 
[^illustrates both the lattice geometries implemented in ABAQUS and the experiment. 

4.1 ABAQUS simulation 

We have compared the wave field for three cases: the first for a homogeneous lattice in the 
absence of any void, the second in the presence of a void, and the third in the presence of a void 
surrounded by our specially designed invisibility cloak. The simulations have been performed 
using a parametric python script for ABAQUS, run by means of MATLAB. We have computed 
the steady-state frequency response of each lattice using the Dynamic/Explicit package already 
implemented in ABAQUS. Since the elastic lattices are constructed from thin elastic ligaments, 
we have performed the simulations employing 3D beam elements. 

The cloak design is an approximate one, in the sense that it works well for a certain range of 
frequencies. Using the simulations we were be able to determine the range of frequencies over 
which the cloak is effective. An out-of-plane vibration in the range between 100 Hz and 250 Hz 
has shown the predicted cloaking action. In particular, in Fig. [T]we show the cases of 120 and 
240Hz, where the presence of the cloak leads the holed plate to behave in a similar way to an 
homogenous lattice without a void. 

4.2 The experimental implementation 

The model, designed theoretically and implemented in ABAQUS, has also been verified experi¬ 
mentally (the geometry is as in Fig. [^. The experiment has been performed at the ‘Instabilities 
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Fig. 6: The comparison between the geometry of the real cloaked lattice drilled by a milling machine (A) and a 
frontal view (B) and a 3D view (C) of the geometry implemented in ABAQUS. The zoomed area of the part (A) 
shows the physical realisation of the square cloak embedded in a square lattice. By observing the enlarged detail 
of the square cloak, the different cross sectional areas of the ligaments are clearly visible and chosen to match 
the required Di and D 2 principal flexural rigidities. 


Lab’ of the University of Trento. In the current setup, we use the Hooke-Chladni-Faraday vi¬ 
sualisation described earlier in Section [2l 

Three structured flexural systems, namely a homogenous lattice, a lattice with a hole, and 
a lattice with a cloaked hole, have been produced by drilling polycarbonate plates (white 2099 
Makrolon UV from Bayer) with an EGX-600 Engraving Machine (accuracy 0.01mm, by Roland). 
The mechanical properties of the polycarbonate, namely elastic modulus (E), density (p) and 
Poisson’s ratio (z/), are respectively 2350 MPa, 1200 kg/m^ and 0.35. The lattice cloak as well 
as the geometrical and material parameters are chosen consistently, according to the same data 
set stored in the SOLIDWORKS file, as illustrated in Fig. Ih, which shows the experimental 
sample. Thin transparent film was used to cover the flexural lattice system to enable the use 
of powder for the Hooke-Chladni-Faraday visualisation. 

The lattices, externally measuring 600mm by 400mm, have ligaments of constant cross- 
section (1.75 mm by 2.5mm) outside the cloak, whereas within the cloak the ligaments have 
variable width and height according to the analytical algorithm of Section 

The lattices, constrained by clamps on the two shorter sides and having the other two free, 
have been excited by using a TIRA Vibrations Test System TV51144 and BOSE ElectroForce 
3300 Series II, connected to the left clamp of the lattices, as shown in the experimental arrange¬ 
ment Fig. 

The maximum amplitude (1 mm peak to peak) of the sinusoidal displacement has been 
imposed by the oscillating clamp connecting the Vibrations Test System to an NI CompactRIO 
system, interfaced with Lab VIEW 2014 (National Instruments). 
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Fig. 7: The vibration apparatus employed in the experiments (A) the details of the constraints on both sides of 
the plate, namely a shaking clamp (B) and a rigid clamped (C). 


4.3 The experimental results and comparison with the computational model 

The qualitative assessment of the effect of the cloak has been carried out using the Hooke- 
Chladni-Faraday technique that shows the positions of the nodal lines of the vibrating plate. 

The boundary conditions are chosen consistently, both in the numerical simulation, and the 
physical experiment: one side of the plate is rigidly clamped, while the opposite side is attached 
to a moving clamp and excited by applying a time-harmonic displacement. The remaining two 
sides are traction free, i.e. the bending moments and the transverse forces are equal to zero at 
these boundaries. 

Photos have been taken with a Sony NEX 5N digital camera (equipped with 3.5-5.6/18-55 
lens, optical steady shot manufactured by Sony Corporation) and with a Nikon D200 digital 
camera (equipped with a AF-S micro Nikkor lens 105 mm 1:2.8G ED). 

Fig. shows the comparison between the experiment and the numerical simulations in the 
case of an applied displacement with a frequency of 120 Hz. 

In particular, in the upper part of Fig. [^we see the unperturbed plate without a hole, with 
the one-dimensional frequency response pattern shown. In the middle part of Fig. [^we show 
that plate containing an uncloaked square hole, with a traction free boundary. 

We observe similar shadow regions both, in the ABAQUS simulation, and the physical 
experiment with a real structured plate. 

Finally, we demonstrate in the lower part of Fig. [^that the shadow region is significantly 
reduced when the cloaking region, depicted in Fig. [^, is introduced around the hole, so that 
the wave pattern outside the hole becomes almost fiat, as expected. It is clear that there is 
good qualitative agreement between the experimental and numerical results. In fact, in the 
case of a cloaked void, the nodal lines showing the incident field represented by a plane wave 
are almost straight, similar to the case of the homogeneous lattice. On the other hand, in the 
absence of the cloak, the pattern lines appear deeply infiuenced by the hole. In this case we 
observe rounded nodal lines that differ significantly from the straight wavefronts observed for 
the homogenous lattice. 

5 Concluding remarks 

In this paper, we have presented a proof of concept design for a square invisibility cloak. Having 
constructed the cloak, we proceeded to examined its effectiveness using both, in a computational 
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ABAQUS model, as well as in real physical experiment. This novel design, proposed in an 
earlier theoretical paper m, was appealing due its simplicity and elegance, which made an 
experimental implementation feasible. The regularisation introduced into this design of cloak 
enables careful and precise implementation of boundary conditions on the interior boundary of 
the cloak in both the numerical simulations and the experiments. 

The approximate cloak presented here, proves to be efficient within a predicted frequency 
range, but the results become frequency sensitive as the frequency of the incident wave increases. 
This effect has been expected, and similar phenomena of the high frequency sensitivity were 
noted in m- 

The range of applications of the proposed cloaking device is wide and it covers, in partic¬ 
ular, earthquake resistant systems, as well as novel designs of foundations of civil engineering 
structures. 

Acknowledgments The authors would like to thank Prof Davide Bigoni for stimulating dis¬ 
cussions and valuable suggestions on the text of the paper, and for his strong support. D.M. 
and A.B.M. gratefully acknowledge financial support from the ERC Advanced Grant insta¬ 
bilities and nonlocal multiscale modelling of materials’ FP7-PEOPLE-IDEAS-ERC-2013-AdG 
(2014-2019) and provision of experimental facilities via the instabilities Lab’ of the Univer¬ 
sity of Trento (http://www.ing.unitn.it/dims/ssmg/). N.V.M gratefully acknowledge financial 
support from European FP7 - INTERCER-2 project (PIAP-GA-2011-286110-INTERCER2). 
D.J.C. thanks the EPSRC for support through research grant EP/J009636/1. 

References 

[1] D. Schurig, J.J. Mock, B.J. Justice, S.A. Cummer, J.B. Pendry, A.F. Starr, and D.R. 
Smith (2006). ‘Metamaterial electromagnetic cloak at microwave frequencies’. Science, 314, 
977980. (doi:10.1126/ science. 1133628) 

[2] R. Fleury and A. A1 (2014). ‘Cloaking and invisibility: A review’. Forum for Electromag¬ 
netic Research Method and Applications Technologies (FERMAT), vol.l, no. 9. 

[3] G. W. Milton, M. Briane, and J. R. Willis (2006). ‘On cloaking for elasticity and physical 
equations with a transformation invariant form’. New J. Phys., 8 , 248 

[4] G.W. Milton, J.R. Willis (2007). ‘On modifications of Newton’s second law and linear 
continuum elastodynamics’, Proc. R. Soc. A, 463, 855-880. (doi:10.1098/rspa.2006.1795) 

[5] A.N. Norris, A.L. Shuvalov (2011). ‘Elastic cloaking theory’. Wave Motion, 48, 525-538 

[6] M. Brun, S. Guenneau, and A.B. Movchan (2009). ‘Achieving control of in-plane elastic 
waves’App/. Phys. Lett., 94, 061903. 

[7] M. Brun, D.J. Colquitt, I.S. Jones, A.B. Movchan, and N.V. Movchan. ‘Transformation 
cloaking and radial approximations for fiexural waves in elastic plates’. New J. Phys., 16, 
093020 (doi:10.1088/1367-2630/16/9/093020) 

[8] M. Farhat, S. Guenneau, S. Enoch, and A. B. Movchan (2009). ‘Cloaking bending waves 
propagating in thin elastic plates’, Phys. Rev. B, 79, 033102 

[9] I.S. Jones, M. Brun, N.V. Movchan, and A.B. Movchan (2015). ‘Singular perturbations 
and cloaking illusions for elastic waves in membranes and Kirchhoff plates’, Int. J. Solids 
Structures, 69-70, 498-506 

[10] D.J. Colquitt, I. S. Jones, N.V. Movchan, A.B. Movchan, M. Brun, and R.C. McPhedran 
(2013). ‘Making waves round a structured cloak: lattices, negative refraction and fringes’, 
Proc. R. Soc. A, 469, 20130218 


11 


[11] D.J. Colquitt, M. Brun, M. Gei, A.B. Movchan, N.V. Movchan, and I.S. Jones (2014). 
‘Transformation elastodynamics and cloaking for flexural waves’, J. Mech. Phys. Solids^ 
72, 131-143. 

[12] S. Guenneau, R.C. McPhedran, S. Enoch, A.B. Movchan, M. Farhat, N.P. Nicorovici 
(2010). ‘The colours of cloaks’, J. Opt, 13, 024014. (doi:10.1088/2040-8978/13/2/024014) 

[13] Parnell W. J. 2012 ‘Nonlinear pre-stress for cloaking from antiplane elastic waves’, Proc. 
R. Soc. A, 468, 563-580 

[14] W.J. Parnell, A.N. Norris, and T. Shearer (2012). ‘Employing pre-stress to generate finite 
cloaks for antiplane elastic waves’, Appl. Phys. Lett., 100 , 171907. 

[15] A. N. Norris and W. J. Parnell (2012). ‘Hyperelastic cloaking theory: Transformation 
elasticity with pre-stressed solids’, Proc. R. Soc. A, 468, 2146 

[16] A.C. Hladky-Hennion, J.O. Vasseur, G. Haw, C. Croenne, L. Haumesser, and A. N. Norris 
(2013). ‘Negative refraction of acoustic waves using a foam-like metallic structure’, Appl. 
Phys. Lett., 102 , 144103 

[17] N. Stenger, M. Wilhelm, and M. Wegener (2012). ‘Experiments on Elastic Cloaking in Thin 
Plates’, Phys. Rev. Letters, 108, 014301. 

[18] T. Biickmann, M. Thiel, M. Kadic, R. Schittny, and M. Wegener (2014). ‘An elasto- 
mechanical unfealability cloak made of pentamode metamaterials’, Nat Commun. 5:4130. 
(doi: 10.1038/ncomms5130) 

[19] T. Biickmann, M. Kadic, R. Schittny, and M. Wegener (2015). ‘Mechanical cloak design by 
direct lattice transformation’, PNAS, 112 (16), 4930-4934. (doi: 10.1073/pnas.1501240112) 

[20] R.V. Craster, J. Kaplunov, and J. Postnova (2010). ‘High-Frequency Asymptotics, Ho¬ 
mogenisation and Localisation for Lattices’, Quarterly Journal of Mechanics and Applied 
Mathematics, 63, 497-519. 

[21] D.J. Colquitt, R.V. Craster, and M. Makwana (2015). ‘High frequency homogenisation for 
elastic lattices’. Quarterly Journal of Mechanics and Applied Mathematics, 68, 203-230. 

[22] A. Colombia, P. Roux, S. Guenneau, and M. Rupin (2015). Directional cloaking of flexural 
waves in a plate with a locally resonant metamaterial J. Acoust. Soc. Am. 137, 1783-1789. 

[23] N.V. Movchan, R.C. McPhedran and A. B. Movchan (2011). ‘Flexural waves in structured 
elastic plates: Mindlin versus bi-harmonic models’, Proc. R. Soc. A, 467, 869-880. 

[24] S. G. Haslinger, N.V. Movchan, A. B. Movchan, and R.C. McPhedran (2012). ‘Transmis¬ 
sion, trapping and filtering of waves in periodically constrained elastic plates’, Proc. R. 
Soc. A, 468, 76-93. 

[25] D. Evans, R. Porter (2007). ‘Penetration of flexural waves through a periodically con¬ 
strained thin elastic plate in vacuo and floating on water’. Journal of Engineering Mathe¬ 
matics, 58 (1-4), 317-337. 

[26] P. Chen, Z. Luo, S. Giiven, S. Tasoglu, A.V. Ganesan, A. Weng, and U. Demirci, (2014). 
‘Microscale Assembly Directed by Liquid-Based Template’, Adv. Mater. 2014, 26, 5936- 
5941. 

[27] Y. Ochiai, T. Hoshi, and J. Rekimoto (2014). ‘Three-Dimensional Mid-Air 

Acoustic Manipulation by Ultrasonic Phased Arrays’, PLoS ONE 9(5): e97590 

doi: 10.1371/journal.pone.0097590. 

[28] E.F.F. Chladni (1787). ‘Entdeckungen iiber die Theorie des Klanges’, Leipzig: Weidmanns 
Erben und Reich. 

[29] L. Rayleigh (1878) ‘The Theory of Sound’, vol. 1 & 2, The Macmillan Company 1877, 
London. 


12 



[30] L. Rayleigh (1911). ‘On the calculation of Chladni’s figures for a square plate’, Philosophical 
Magazine Sixth Series 22 225-229. 

[31] W. Ritz (1908). ‘Uber cine neue Methode zur Losung gewisser Variationsprobleuie der 
mathematischen Physik’, Journal fiir Reine und Angewandte Mathematik 135 1-61 

[32] W. Ritz (1909). ‘Theorie der Transversalschwingungen einer quadratische Platte uiit freien 
Randern’, Annalen der Physik 28 737-786. 

[33] J.B. Pendry, D. Schurig, D.R. Smith (2006). ‘Controlling Electromagnetic Fields’, Science 
312 (5781): 1780-1782. 


13 



